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Abstract. Taking into account the electromagnetic and gravitational form factors, calculated from a new 
set of t-dependent GPDs, a new model is built. The real part of the hadronic amplitude is determined only 
through complex s. In the framework of this model the quantitative description of all existing experimental 
data at 52.8 < y/s < 1960 GeV, including the Coulomb range and large momentum transfers (0.0008 < \t\ < 
9.75 GeV 2 ), is obtained with only 3 fitting high energy parameters. The comparison with the preliminary 
data of the TOTEM Collaboration at an energy of 7 TeV is made. 

PACS. 13.40.Gp , 14.20.Dh, 12.38.Lg 



1 Introduction 

TlJeLdynamics f strong interactions finds its most com- 
plete representation in elastic scattering at small angles. 
Only in this region of interactions we can measure the 
basic properties of the non-perturbative strong interac- 
tterHwhich defines the hadron structure: the total cross 
section, the slope of the diffraction peak and the param- 
eter p(s,i). Their values are connected, on the one hand, 
with the large-scale structure of hadrons and, on the other 
hjjnd, with the first principles which lead to the theorems 
OT^jAie behavior of the scattering amplitudes at asymptotic 
energies [TJ[2]. 

J 'herp are indeed many different models for the de- 
s^ption of hadron elastic scattering at small angles [3l 
HJ^Nphey lead to the different predictions for the struc- 
ture of the scattering amplitude at asymptotic energies, 
w&eje the diffraction processes can display complicated 
features [5]. This concerns especially the asymptotic uni- 
tarity bound connected with the Black Disk Limit (BDL) 
,[6feJt was assumed in the Chow- Yang model [3[8] that the 
hadron interaction is proportional to the overlapping of 
the matter distribution of the hadrons and in the Wu and 
Yang [7 that the matter distribution is proportional to 
the charge distribution of the hadron. Then many models 
used the electromagnetic form factors of hadron, but, in 
most part they change its form to describe the experimen- 
tal data, as was made in the famous Bourrely-Soffer-Wu 
model [9] . The parameters of the obtained form factor are 
determined by the fit of the differential cross sections. 

Now we present a model, based on the assumption that 
the hadron interaction is sensitive to the generalized par- 
ton distributions (GPDs) whose moments can be repre- 
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sents over momentum transfer in the form of two differ- 
ent distributions: charge and matter, separately. Hence, 
this model uses the exact electromagnetic and gravitation 
form-factors determined by one function - generalized par- 
ton distributions (GPDs). So both form factors are inde- 
pendent of the fitting procedure of the differential cross 
sections. Note that the form of the GPDs is determined, 
on the one hand, by the deep-inelastic processes and, on 
the other hand, by the measure of the electromagnetic 
form factor from the electron-nucleon elastic scattering. 
Hence, the form of the electromagnetic form factor (first 
momentum of GPDs) determines the form of the second 
form factor (second momentum of GPDs). This picture 
is supported by the good description of the experimen- 
tal data in the Coulomb-hadron interference region and 
large momentum transfer at high energies by one ampli- 
tude with a few free parameters. The impact of the hard 
pomeron contribution to the elastic differential cross sec- 
tions is very important for understanding of the properties 
of the QCD in the non-perturbative regime [TU]. Note that 
the real part of the hard pomeron is essentially large than 
the real part of the soft pomeron. Now in [TT] it is sug- 
gested that such a contribution can be explained by the 
preliminary result of the TOTEM Collaboration [T^] on 
the elastic proton-proton differential cross sections. In our 
model, the real part of the hadronic amplitude is deter- 
mined only through complex s satisfying the cross sym- 
metric relation. In the framework of this model, the quan- 
titative description of all existing experimental data at 
52.8 < y/s < 1960 GeV, including the Coulomb range and 
large momentum transfers 0.0008 < \t\ < 9.75 GeV 2 , is 
obtained with only 3 fitting high energy parameters. The 
comparison of the predictions of the model at 7 TeV and 
preliminary data of the TOTEM collaboration are shown 
to coincide well. There is some small place, especially in 
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the region of the diffraction dip, for the small correction 
contributions which are determined by the odderon, and 
possibly the spin-dependent part of the scattering ampli- 
tude which gives a small contribution at large momen- 
tum transfer. In the framework of the model, only the 
Born term of the scattering amplitude is introduced. Then 
the whole scattering amplitude is obtained as a result of 
the unitarization procedure of the hadron Born term that 
is then summed with the Coulomb term. The Coulomb- 
hadron interference phase is also taken into account. The 
essential moment of the model is that both parts of the 
Born term of the scattering amplitude have the positive 
sign, and the diffraction structure is determined by the 
unitarization procedure. 

The electromagnetic and hadronic parts of the elastic 
scattering amplitudes used in the model are presented in 
the second and third sections. In the fourth section, we 
introduce the hadron form factors obtained from the first 
and second momenta of GPds. Our fitting procedure and 
the description of the high energy differential cross sec- 
tions of the proton-proton and proton-antiproton scatter- 
ing are presented in the fifth and sixth sections. Also, we 
stretch our model on 7 TeV and compare the model calcu- 
lations with the preliminary data of the TOTEM Collab- 
oration. Then, the model calculations are compared with 
the experimental data at low energies y/s — 30 GeV. In 
the seventh section and in the conclusion, we show and 
discuss the model calculations for the total cross sections 
and the value of p(s,t) obtained in the framework of the 
model. 



2 Electromagnetic part of the hadron 
scattering amplitude 

The differential cross sections of nucleon-nucleon elastic 
scattering can be written as the sum of different helicity 
amplitudes: 



where the form factors are: 
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The total helicity amplitudes can be written as ^(s, t) — 
F^{s,t)+Ff m (s,t)e^ s '^ , where F t h (s,t) comes from the 
strong interactions, F° m (s, t) from the electromagnetic in- 
teractions and ip(s,t) is the interference phase factor be- 
tween the electromagnetic and strong interactions [T31H41 
[T5] , For the hadron part the amplitude with spin- flip is 
neglected in this approximation, as usual at high energy. 

The electromagnetic amplitude can be calculated in 
the framework of QED. In the high energy approximation, 
it can be obtained |T7j for the spin-non-flip amplitudes: 
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with k relative to the anomalous magnetic moment, and 
Gd(t) has the conventual dipole form 

Gd(t) = 1/(1 -t/0.71) 2 . (5) 



3 Main hadronic amplitude 

The model is based on the representation that at high en- 
ergies a hadron interaction in the non-perturbative regime 
is determined by the reggenized-gluon exchange. The cross- 
even part of this amplitude can have 2 non-perturbative 
parts, possible standard pomeron - Pznp) and cross-even 
part of the 3-non-perturbative gluons (P3« p ). The interac- 
tion of these two objects is proportional to two different 
form factors of the hadron. This is the main assumption 
of the model. Of course, we can not insist on the origin 
of the second term of the scattering amplitude. It can be 
a different nature. However, in any case, it has the cross- 
even properties and positive sign. The second important 
assumption is that we chose the slope of the second term 4 
time smaller a slope of the first term, by analogy with the 
two pomeron cut. Both terms have the same intercept. 

The form factors are determined by the General parton 
distributions of the hadron (GPDs). The first form factor, 
corresponding to the first momentum of GPDs is the stan- 
dard electromagnetic form factor - G(t). The second form 
factor is determined by the second momentum of GPDs 
-A{t). The parameters and i-dependence of the GPDs are 
determined by the standard parton distribution functions, 
so by the experimental data on the deep inelastic scatter- 
ing, and by the experimental data for the electromagnetic 
form factors (see [T5]). 

Hence, the Born term of the elastic hadron amplitude 
can be written as 

F*° rn (s,t) = h x G 2 it) F a (s,t) (l + n/s - 5 ) (6) 
+ h 2 A 2 it) F b (s,t) (l + r 2 /J - 5 ) 

where F a (s,t) and Ff,(s,t) has the standard Regge form 

F a (s,t) =s ei e B ^ *; F b (s, t) = J ei e B ^^ *, (7) 

with G(t) — Gsit) is the Sachs electric form factor rela- 
tive to the first moment of GPDs and A{t) relative to the 
second moment of GPDs. 

Li 4ml - (1 + k) t 
GH) = TTT^ " , ' (8) 
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Fig. 1. The hadron form factor A(t) - (the second momentum 
of GPDs) in the t representation (hard line - L 2 — 2.0 GeV 2 ; 
dashed line - L 2 = 1.8 GeV 2 ; circles - the numerical calculation 
of the integral over a:) 



with the parameters: L\ = 0.71 GeV 2 ; L\ = 2 GeV 2 . 

7«o; so - 1 GeV 2 . (10) 



s = se -"/ 2 / 



The slope of the scattering amplitude has the standard 
logarithmic dependence on the energy. 



B(s) = a' ln{s). 



(11) 



this i-dependence (TB] based on the original Gaussian form 
corresponding to that of the wave function of the hadron. 
It satisfies the conditions of non-factorization, introduced 
by Radyushkin, and the Burkhardt condition on the power 
of (1 — x) n in the exponential form of the ^-dependence. 
With this simple form we obtained a good description of 
the proton electromagnetic Sachs form factors. Using the 
isotopic invariance we obtained good descriptions of the 
neutron Sachs form factors without changing any param- 
eters [16] . 

We shall use this model of GPDs to obtain the sec- 
ond momentum form factor of the nucleon. Taking instead 
of the electromagnetic current the energy-momentum 
tensor T^ v together with a model of quark GPDs, one can 
obtain the gravitational form factor of fermions 

J ^ dx x [H(x, A 2 ,0± E(x, A 2 ,0] (16) 
= A q {A 2 )±B q (A 2 ). 

For 5 = one has 

/ dx x[H(x,t)±8(x,t)] =A q (t)±B q (t). (17) 
Jo 

The integration of the second momentum of GPDs over 
x gave the momentum-transfer representation of the form 
factor (see Fig.l). We approximate this by the dipole form 



with a' = 0.24GeV~ 2 . 

The final elastic hadron scattering amplitude is ob- 
tained after unitarization of the Born term. So in first, we 
have to calculate the eikonal phase 



A(t) = Li/(L 2 2 -t) 2 . 
with the parameter L\ = 2.0 GeV 2 . 



(18) 



x(s,b) = J d 2 q e *-« F**» ( s , q 2 ) , (12) 5 Fitting procedure 

and then obtain the final hadron scattering amplitude 

F h (s,t) =is [ b J (bq) r(s,b) db. (13) 



r(s,b) = l-cxp[~ X (s,b)} 



(14) 



All these calculations are carried out by the FORTRAN 
program. 

4 Hadron form factors 

As was mentioned above, all the form factors are obtained 
from the GPDs of the nucleon [IB] . The electromagnetic 
form factors can be represented as first moments of GPDs 

F?{t)= f dxH q (x,t); F|(t) = / dx £ q (x,t), (15) 



The model has only 3 high energy fitting parameters and 
2 low energy parameters, which reflect some small contri- 
bution coming from the different low energy terms, (see 
Table 1). We take all existing experimental data in the 
energy range 52.8 < y/s < 1960 GeV and the region of 
the momentum transfer 0.0008 < -t < 9.75 GeV 2 of 
the elastic differential cross sections of proton-proton and 
proton- antiproton data [50,21 . So we include the whole 
Coulomb-hadron interference region where the experimen- 
tal errors are remarkably small. We do not include the data 
on total cross sections atot{s) and p(s), as their values were 
obtained from the differential cross sections especially in 
the Coulomb-hadron interference region. Including such 
data decreases % 2 . We also do not include the interpolated 
and extrapolated data of Amaldi |22j . 

In the fitting procedure we calculate the minimum in 
Si=i X? related with the statistical errors of. The sys- 
tematical errors are taken into account by the additional 
normalization coefficient ilk for the k series (the experi- 
ment) of the experimental data 



following from the sum rules [181119) . 

Recently, there were many different proposals for the 
t dependence of GPDs. We introduced a simple form for 



i=l 



n k Ej(s,t) - Tj(s,t) 



(19) 
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high energy parameters low energy parameters 



hi, GeV" 2 


h 2 , GeV" 2 


ei 


n, GeV 


r 2 , GeV 


0.23 ±0.01 


0.74 ±0.01 


0.11 ±0.01 


11.95 ±0.5 


-5.9 ±0.8 



Table 1. The basic parameters of the model are determined by fitting experimental data. 




Fig. 2. da/dt at y/s = 52.8 GeV and at small t for pp (left) and for pp (right). 



where Ti(s, t) are the theory predictions and nkE{(s, t) are 
the data points allowed to shift by the systematical error 
of the fc-experiment (see, for example p?3"l[2~4"] . 

In the region of the small momentum transfer the sys- 
tematic errors are of an order of 2% -j- 5%. For most part 
the additional normalization are in the region 0.95 ± 1.05. 
At large momentum transfer the order of the systematical 
errors is 10% ± 20%. In this case, the additional normal- 
ization is situated in the region 0.8 ± 1.2. 

For the non- normalized experimental data of the U AA/2 
Collaboration [25] , which have very small statistical er- 
rors, we take the normalization determined in 26 . Our 
correction normalization is obtained from the fitting pro- 
cedure in this case nuA43 = 0.95. 

As a result, one obtains \\ /N ~ 1.8 where N = 975 
is the number of experimental points. Of course, if one 
sums the systematic and statistical errors, the X 2 /N de- 
creases, to 1.4. Note that the parameters of the model are 
energy-independent. The energy dependence of the scat- 
tering amplitude is determined only by the single intercept 
and the logarithmic dependence on s of the slope. 

Note that there are some separate points ( n = 17) at 
the different energies and momentum transfer which give 
Y^n=i Xn = 260. However, we do not remove such points. 

6 Description of the differential cross sections 

The differential cross sections for proton-proton elastic 
scattering at y/s = 52.8 GeV are presented in Fig. 2 (left 
panel) and 3 (top panel). At this energy there are exper- 
imental data at small (beginning at —t = 0.001 GeV 2 ) 
and large (up to —t = 10 GeV 2 ) momentum transfers. 
The model reproduces both regions and provides a quali- 
tative description of the dip region at — t ~ 1.4 GeV 2 , for 



yfi = 53 GeV 2 and for y/s = 62.1 GeV 2 (Fig.3(top) and 
Fig.4(lctt panels)). 

Now let us examine the proton-antiproton differen- 
tial cross sections (see Fig.2(right panel)). In this case 
at small momentum transfer the Coulomb-hadron inter- 
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Fig. 3. da/dt at y/s = 52.8 GeV at large \t\ for pp (top) and 
for pp(bottom panel) 
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Fig. 4. da/dt at y/s = 62.1 GeV at lar f 

ference term plays an important role and has the oppo- 
site sign. The model describes the experimental data well. 
Slightly worse than pp is the description of pp of differ- 
ential cross sections at y/s = 53 GeV, especially in the 
diffraction minimum (Fig.3(bottom panel)). 

Maybe, this shows an additional odderon contribution. 
Note that at y/s = 62.2 GeV for pp scattering the descrip- 
tion of the differential cross sections is essentially better 
(Fig. 4). 




10-5 I I 

0.0 0.5 1.0 1.5 2.0 

-t (GeV2) 

;e \t\ for pp (left) and for pp(right panel) 

In Fig. 5, the description of the proton-antiproton scat- 
tering at y/s — 541 GeV and at y/s = 1800 GeV is shown. 
In this cases, the Coulomb-hadron interference term is 
large, especially at y/s = 541 GeV as t is very small. The 
good description of the experimental data shows that the 
energy dependence of the real part of scattering amplitude 
obtained in the model corresponds to the real physical sit- 
uation. 





Fig. 5. da/dt for pp elastic scattering at small momentum Fig. 6. da/dt for pp elastic scattering at large t, at y/s = 
transfer, at y/s = 541 GeV (top) and y/s = 1800 GeV (bottom) 546 GeV (top) and y/s = 630 GeV (bottom panel) 
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Figures 6 and 7 show the description of the experimen- 
tal data at larger momentum transfers for yfs = 546 GeV 2 
and yfs — 630 GeV 2 and for Tevatron energies yfs = 1800 
GeV and yfs — 1960 GeV. It is clear that the model leads 
to a good description of these data in the region of the 
diffraction minimum without taking into account the odd- 
eron contribution. Hence, it is shown that very likely the 
intercept of the odderon is near 1. 

On basis of this fit of the experimental data at 52.8 < 
sfs < 1960 GeV and 0.0008 < \t\ < 9.75 GeV 2 we ob- 
tained the fitting parameters (see Table 1). Taking into 
account these values of the parameters we extend the 
scope of the model and calculate the differential cross sec- 
tions at 7 TcV for pp elastic scattering. In Fig. 8, the 
comparison of the model calculations with the parame- 
ters, based at the fit of the existing experimental data at 
52.8 < yfs < 1960 GeV, with the preliminary data of the 
TOTEM Collaboration are shown. Except the size of the 
diffraction minimum the coincidences are remarkable. Of 
course, if we include in the model some different correction 
terms, like odderon, the value of the fitting parameters of 
the model will be slightly change. However, we think that 
the basic properties of the model will not change in future. 

Now let us see how it can be extend the scope of the 
model on a low energy . The calculation of the model of 
da/dt for pp and pp elastic scattering at yfs = 30.6 GeV 
are compared with the experimental data on Fig. 9 and 
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Fig. 8. The comparison of the model calculations with the 
parameters, based on the fit of the experimental data at 52.8 < 
yfs < 1960 GeV, with the preliminary data of the TOTEM 
Collaboration at yfs = 7 TeV. 

Fig. 10. There is a very good description at small momen- 
tum transfer for both the reactions - pp and pp. At large 
t the model reproduces the differential cross sections only 
qualitatively. The position of the diffraction minimum cor- 
responds to the experimental data. However, the form of 
the diffraction minimum obviously requires some addi- 
tional small correction terms, possibly by the odderon con- 





Fig. 7. da/dt for pp elastic scattering at large t, at yfs = 
1800 GeV (top) and yfs = 1960 GeV (bottom panel) 



Fig. 9. The model predictions of da/dt for pp elastic scattering 
at at y/s = 30.6 GeV. 
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Fig. 10. The model predictions of da/dt for pp elastic scattering at at yfs = 30.6 GeV. 



tributions. This situation repeated the problem of describ- 
ing the form of the diffraction minimum at y/s = 53 GeV 
for pp elastic scattering. 
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Fig. 11. p(s,t) for pp elastic scattering at t = (top) (cir- 
cle and triangles are the experimental data of pp and pp elas- 
tic scattering, respectively); and at small momentum transfer 
(bottom) at yfs = 52.8 GeV (tiny-dash line), yfs = 541 GeV 
(dash line), yfs = 7 TeV (hard line), and = 14 TeV 
(long-dash line). 




Fig. 12. atot(s) are calculated in the model (hard line), 
(circle and triangles are the experimental data of pp and pp 
elastic scattering, respectively. 



7 Energy dependence of p(s,t) and a to t{s) 

The ratio of the real part to the imaginary part of the 
elastic scattering hadronic amplitude 

p(s,t) = ReF h (s,t)/ImF h (s,t) 

is very important as it reflect the t-dependence of the both 
parts of the scattering amplitude, which are connected one 
to the other through the integral dispersion relations. The 
validity of this relation can be checked at LHC energies. 
The deviation can point out the existence of a fundamen- 
tal length at TeV energies [271128] , Usually, the value of 
p(s, t = is assumed to be small and to vary little with 
t: p(s,t) w 0.14. The differential cross sections at small 
momentum transfer \t\ < 0.05 GeV 2 , so-called Coulomb- 
hadron interference region, are determined by the inter- 
ference of the Coulomb amplitude with the real part of 
the hadron amplitude. Hence, the s and t dependence of 
the real part of the hadron amplitude, which is reflected 
in p(s,t), will determine the form of the differential cross 
sections. 

In the model, the real part of F/,(s,t) is determined 
only by the complex cross symmetric form of energy s = 



8 



O. V. Selyugin: GPDs of the nucleons and elastic scattering at high energies 



s exp[—iir/2\. No other artificial function or some param- 
eters impact the form or s and t dependence of the real 
part. Despite such simplicity, the model sufficiently well 
describes the experimental data in the Coulomb-hadron 
region of momentum transfer and in a wide energy re- 
gion (see figs. 2,5,9). The calculated p(s,t = 0) in the 
model is shown in Fig. 11 (top panel). In most part it co- 
incides with p obtained by the COMPETE Collaboration 
[2"M5tj] . We can see that the p(s,t) reaches its maximum 
0.135 at approximately 1 TeV and then slowly decreased 
up to 0.127 at y/s = 20 TeV. Note that the model takes 
into account the only cross-symmetric part of the scat- 
tering amplitude. This reflects that at low energy it, in 
most part, coincides with the experimental data of the 
proton-proton scattering. The t dependence of p, shown 
in Fig. 11 (bottom panel), is interesting. We can see that 
after y/s = 1 TeV, where p reaches its maximum, its t- 
dependence is changed. If the maximum of p is small de- 
creasing, his slope will be larger and larger. It is interesting 
to compare this figure with a figure like this in our phe- 
nomenological analysis of p at high energies [31] . The size 
of p is essentially larger, but the dynamic change of p with 
s and t is practically the same. It is related with the sat- 
uration regime which impacts on the s and t dependence 
of p(s,t) at high energies. 

As the data on the total cross section cr tot (s) and p(s, t) 
were not included in our fitting procedure, let us now cal- 
culate these values in the framework of the model. The 
energy dependence of the atot(s) is shown in Fig. 12. Our 
model calculations coincide with the existing experimen- 
tal values sufficiently well. The model calculations gives 
a t ot = 43.46, 61.7, 75.76 mb for the energies - y/s = 
50, 540, 1800 GeV. 

8 Conclusion 

We present the new model of the hadron-hadron inter- 
action at high energies. The model is very simple from 
the viewpoint of the number of parameters and functions. 
There are no any artificial functions or some cuts which 
bound the separate parts of the amplitude by some re- 
gion of momentum transfer. The model describes all high 
energy data sufficiently well, except the form of the diffrac- 
tion minimum of the proton-antiproton scattering at y/s = 
53 GeV. As we know, it is the only model which de- 
scribes all available high energy data in the Coulomb- 
hadron region and large momentum transfer. The energy 
dependence of the differential cross sections is determined 
by only one intercept A = 1 + e with e = 0.11 and 
the a' = 0.24 ln[s\. One of the most remarkable prop- 
erties is that the real part of the hadron scattering ampli- 
tude is determined only by complex energy s that satisfies 
the crossing-symmetries. The differential cross sections at 
small momentum transfer are determined, in most part, by 
the Coulomb-hadron interference term. This is essentially 
caused by the s and t dependence of the real part of the 
hadronic amplitude. A good description of the differential 
cross section (see figs 2, 5, 9) confirms the determination 
of the real part in the model. 



However, the most important advantage of the model 
is that it is built on some physical basis - two form fac- 
tors which are calculated from GPDs. The behavior of the 
differential cross section at small t is determined, in most 
part, by the electromagnetic form factors; and at large t, 
by the matter distribution (calculated in the model from 
the second momentum of the GPDs), as was supposed by 
H. Miettinen a long time ago [32] and then used in the the 
work of S. Sanielevici and P. Valin [33 . Of course, we un- 
derstand that this advantage of the model, on the other 
hand, has some disadvantage. Now we have no rigorous 
proof of such a physical picture. However, the best work 
of the model maintains the hope that such a representa- 
tion has a right to exist. 

The model predicts a to t — 95 mb at y/s = 7 TeV. 
The preliminary data of the TOTEM show a to t = 98.3 ± 
0.2 stat - ± 2.8 svst ^ mb [J]. Of course, some correction terms 
(corresponding to the odderon or spin-flip amplitudes) 
have to be included in our model. Now the model does 
not show a contribution of the hard pomeron in the exam- 
ined energy region. The presence of the hard pomeron has 
to give an additional contribution in the real part of the 
scattering amplitude and change the size and form of the 
Coulomb-hadron interference term |31) . Hence, we hope 
that all these additional terms will be determined after 
fitting with the new data of the proton-proton scattering 
at LHC energies. 
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